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SUMMARY 



A mathematical procedure is herein developed tor 
obtaining exact solutions of shear— lag problems in flat 
panels and "box "beams'":' the method is "based on the assump- 
tion that the amount of stretching of the sheets - in" the 
direction perpendicular to the dir ect ion" of" 'e¥¥e5tfal~ 
normal stresses ia negligible. "Explicit solutions, in- 
cluding the treatment of cut— outs, are given for several 
cases and numerical results are' pr esehjbed in grapfcic'^and 
tabular form. The general 'theory is -pre'se'nted in a fori 
from which further solutions can .be" readily obtained. , 
The extens ion "of ' the theory to cover certain c ases"~of non~" 
uniform cross section/is indicated. Although the^s"61u— 
tions are obtained in terms of infinite series, the pres- 
ent developments differ from those previously given in 
that, in practical cases, the series usually converge 90 " ' 
rapidly that sufficient accuracy is afforded by a small 
number of terms. Comparisons are made in several cases" 
"between the present results and the corresponding solu- 
tions- obtained' by approximate procedures devised by 
Reissner and by Euhn and Chiarito. - ■- t-~ 



INTRODUCTION 



The theory given in this paper is to be considered 
as a refinement of approximate methods devised by "R eis is her ~ 
in reference 1 and by Kuhn and ChiaritoTn reference 2 for 
the analysis of shear— lag problems and may thus serve as a 
check on the accuracy of those methods. ""Such problems 
concern the effect of shear deformation oh the state of 
stress in th\n sheets which are analyzed in' the elementary 
theory of the strength of materials without regard to shear 
def ormability . Examples with which this" paper Is concerned 
are: (l) the introduction of concentrated forces into flat 
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sheets, by means of stiffeners, where the forces act in the 
plane of the sheet, and (2) the distribution of stress in 
the cover Bheets of rectangular box beams subjected to 
bending loads: that is, in particular, the determination 
of the effective width of suoh cover sheets. 

One of the assumptions made in the above-mentioned , 
theories, which is retained in the following developments, 
is that stretching of the sheets in the direction perpen- 
dicular to the direction of essential normal stresses can 
be neglected. This assumption, which is made plausible 
by elastic energy considerations (reference 1), is essen- 
tial for the application of the mathematical methods em- 
ployed in most of this paper.. 

No further simplifying assumption, however, is made 
here and the problems considered- are then solved in an 
exact ffianner, as boundary- value problems in the theory of 
plane stress. This work is in contrast with that of 
Reissner, who obtains approximate solutions in substan- 
tially simpler form of a more extensive class of problems 
than the problems considered in the present analysis by 
an application of the principle of least work, and with 
that of Zuhn and Ohiarito, who further simplify the prob- 
lems by Incorporating the resistivity of the sheets in the 
direction of th-e-essential normal stresses into the ef- 
fects of stiffeners in this direction by increasing the 
actual areas of the stiffeners by an amount representing 
an estimated effective sheet width. Kuhn and Chiarito 
then have to deal with sheets that are rigid in the trans- 
verse direction and offer no resistance in the longitu- 
dinal direction but possess finite shear rigidity. 

The following investigation will include the analysis 
of beams with and without cut-outs, and, in this connec- 
tion, it is remarked that while the way was known in which 
series solutions could be obtained for certain problems 
concerning beams without cut-outs, in terms of assumed 
trigonometric developments -in the spanwise direction (ref- 
erences 3 and 4), a~new type of series development is 
obtained herein which makes possible the treatment of both 
types of problem in an analogous way and, in addition, 
possesses the advantage of being much more -adaptable to 
computation than the known type of series solution for tho 
problem of the beam .wi thout cut-out. 

This investigation, conducted at the Massachusetts 
Institute of Technology, was sponsored by, and conducted 
with, financial assistance from the National Advisory 
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CommittpR for Aeronautics. Thanks are due to Professor 
7. Ssissnpr of the Massachusetts Institute of Technology 
for valuable suggestions and criticisms. 

SYMBOLS 

I length of panel or "beam 

w . one— half width of panel or "beam 

x spanwise distancje ^^-l = 

A cr os s— BPCt ional area 

y transverse distanoe 

E externally applied concentrated force 

t sheet thickness 

h height of box "beam 

I y principal moment of inertia of cross section "f two 
side wots , including corner flanges .. _"' 

a normal Stress 

T shear stress 

e normal strain 

«. displacement component in x direction 

v displacement component in y direotion 

E Young's modulus (E = E x ) 

V PoiB son's ratio . . 

Y shear strain * 

G shear modulus ; , 

H - stress function 

c constant 
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H* modified stress function £ H + ffv(x)] 
a ratio of flange and sheet areas 

am = 2tw 




t w. 



n positive integer 

\ auxiliary parameter 

I dimens ionless span coordinate ( j — ~£~J 

t\ d imens i onle s s transverse coordinate 



■ s 




G TTX*\ 



P ratio of moments of inertia of web and sheet about 
transverse principal axis of beam (l w /l g ) 



moment of lnprtia of cross section of two cover sheptB 
about transverse principal beam axis (twh.2) 



p, auxiliary parameter ^ ~ 
p 0 uniformly distributr&d load 

a 0 maximum e le-men-trar y spanwi6e normal stress [ffj, ( I) ] 

6 angle between x— axis and normal to edge of sheet 
Tfi - (2n - I)tt /e w 

L n 2 A 0 

I total principal moment of inertia of beam section 
(I s + I„ + %A ffi h2) 



(yf) 



p particular solution of-^iif f e.r ent ial equatl>n 
k constant 
M moment of applied bending l«f»ads 
P particular solution. of differential equation 
$ stress function 

X,Y separation functions in solut-ion of differential 
equation 

K dimensiinlflss constant occurring in Kuhn theory (with 
subscript R, occurring in He-is s ner ~ trhe or y) 

K auxiliary parameter 
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Subscripts: 

e edge 

m middle 

eff effective 

x spanwise 

y transverse 

n corresponding to n 

b "bending 

a average 



DEFINITION 02 THE PROBLEMS 



The first problem treated concerns the stress dis- 
tribution in a flat panel of length \ and width 2w 
loaded at one end (x = 0) by concentrated axial forces 
introduced into flanges of cross-sectional area JL e 
attached to the longitudinal edges (y = ±v$ and by a 
concentrated axial force introduced into a longitudinal 
stiffener of area A' m along the axis of symmetry (y = 0) 
(See fig. 1,) Suitable conditions, are to be prescribed 
along the edge x = I. For example, (1) the panel may 
be completely clamped along this edge, (2) the flanges 
and the longitudinal may be fixed at x = \ with the 
sheet subject- to displacement, or (3) the flanges may 
be fixed while both the sheet and the stiffener are not 
attached to a support* It will be convenient to consider 
separately loadings that are symmetrical and antisymmet- 
rical with respect to the longitudinal axis of the panel 
(fig. 1) and to obtain solutions for arbitrary loadings 
by superposition. 

The second problem deals with the stress distribu- 
tion in the cover sheets of a d.oubly symmetrical box 
beam x supported in a prescribed manner at one end 
(x = \) and unsupported at the other end (x = 0) and sub- 
jected to a given distribution of bending loads applied 
symmetrically at the sheet-web Junctions (fig. 2). 
There may be flanges at the junctions of the cover 
sheets and side webs and a longitudinal stiffener along 
the center line (y = o) of the cover sheets. The pro- 
cedure will be outlined for a rather general class of 
end-support conditions, while explicit solutions will 
be obtained for the cases when (1) the structure is 
completely clamped along the edge at x = \ (fig. 
2(a)), and (3) the side webs are fixed at x = \, 
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but the cover sheets are not attached to a support (fig. 
2(b)), as would be tho case at the outboard edge of a 
cut-out . 

Finally, the Btress distribution is determined in an 
infinite half-sheet rigid in the direction parallel to 
its straight-edge and loaded by a concentrated force, in 
the plane of the sh«et and normal to the straight-edge, 
introduced into a stiffener attached to the sheet in the 
direction of the applied force (fig. 3(a)). Vhile the 
solution in this case is probably of slight practical 
interest in itself, it was felt that it might serve as 
an indication of the limiting behavior of solutions to 
other related problems. 



BASIC EQUATIONS IN THE THE OB Y OF PLANE STRESS FOE 
AN ORTHOTROPIC MATERIAL RIGID IN ONE DIRECTION 



The following equations mu.gt be satisfied for a 
state of plane st"r~e~s~s~ in an orthotropic medium of uni- 
form thickness; 

(a) The - equilibrium conditions for an element of the 
sheet • • • - • 



3<?x + 3t 



ox 9y 



0 



(la) 



6t 3py 
3x 9y 



0 



(lb) 



(b) The stress- stTain relationships 



cfx = E~ ( ^* ' 



3u 1- 



u xOy> 



(2a) 



y 



W * E7 ( °y - 




(2b) 



3u "9v 1 
5" y + 3x ~ G 



T 



(2o) 



For a medium rigid in the y-direct-ion 
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Ey = CO 



v x = 0 



(3) 



and, from equation (21)), it then follows that 

dv 

v = v(x) 



= 0 



(4) 



If E is written for E x , equations (3a) and (2c) take 
the form " " 



E 



3U 



ou 



G — = T - G- t'(x) 

07 



(5a) 



(5b) 



? rom equation (5) it follows by differentiation that 



If the equilibrium condition of equation (la) ia 
satisfied by means of a stress function H, in terms of 
which 

dH 



Q3C 



(7a) 
(7b) 



it follows from equation (5) that H must satisfy the 
equat ion 

o 5 H G o 5 H 
dxa E dy s 



+ - - = — & T«(x) 



("a) 



The solution of any plane stress problem in such an 
orthotropic medium is thus reduced to the formulation of 
boundary conditions relative to the function H and to 
the solution, of the differential equation (8) subject to 
those boundary conditions. 

With H determined, the stress a y is found from 
eauation (lb) as follows* .. 
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bo. 



y 
by 



St b 3 E 



G- 5 3 H 



bx bx 2 



& v«(x) 



E by 2 
& y v«(x) + .f(x) 



(9) 



while from equation (5) there follows for the displace- 
ment component u, 



E J 



ay 



dx 



(10) 



An interesting .relationship . "between u and H can 
"be established by means of complex variables. If a new 
var iable 



X * = /I X 



is introduced, equation (8) transforms into 



(11) 



a 2 H d 2 n 

+ = — G 



6 Xl 2 by 3 
and equations (5) and (7) give 



d E v 



(13) 



bH _ ^ u 
by bx x 



bH dv rr Su 

+ & = — VE£ — 

ox x dx x _ by 



> 



(13) 



If the function 

K* = H + & v( x) 
is defined, equations (12) and (13) "becom-e 

& 2 H* b 2 H* 



bx x 3 by : 



= '0 



(14) 



(15) 
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ou 

5Xi 



v£g~^ 

oy 



> 



(16) 



Prom equations ( 15 ) and (16), it may "be concluded that 

H* + i u = Idi + iy) (17) 

that is, H + Gv(x) and u are the real and the 

imaginary parts, respectively, of a function. of the com- 



plex variable /— x + iy. 



SUMMARY AND DISCUSSION OF EXPLICIT SOLUTIONS 



In this section are presented explicit, solutions of 
cortain problems, the mathematical derivations of which 
are given in the two sections that follow. Solutions of 
a large class of related problems can be obtained directly 
from the analysis included in those sections. 

Although the results with regard to panels are pre- 
sented in this section only for the limiting case in which 
I = oo , the corr ©spending solutions for panels clamped at 
the end x = I 



■i x. - 

11 w 



by. cosh 



can be obtained in all cases by irepTacing 

/^n 1 

V Gr U W 



cosh 



in the given 



"expressions for the stress functions. 



Unless otherwise specified, the ratio &/E is 
assigned the value 3/8 in the numerical calculations of 
this paper. In case the loading is such that the sheet 
becomes wrinkled, a smaller value of the shear modulus 
must be taken. The nature of the solutions for large val- 
ues of the span— width ratio is such, however, that a 
change in the value of G merely involves magnification 

of the 



x 
w 



scale in the ratio 



y 



G 



G eff ' 



together wltTT a 
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multiplication of the shear Stress values 'by the recipro- 
cal ratio. 



Stresses in a Long Symmetr ically Loaded ITlangod Panel 
without Central Btiffener 
[Panel l] 

A long panel of width 2w is clamped along one end 
(x = oo ) and is loaded at the eni x = 0 by axial forces 
3P acting on flanges of equal cr o e s— s e ct ional ar-ea A 
attached to the edges (y = ± w) of— the sheet (fig. 4(a)). 
The stress function is determined in tire "form 



\ • cos X n sin XnJw ' J T n w I 



H = *< i~T I + 2 > r; ' vn ' M ^ w '- , e v * " w k (18) 

t j^a + 1 w X n (l + a cos 2 X n ) j 

where a, represents the ratio of the flange and shoet area's 
and is given as 

frw — ; 

and the parameters \ n are the ;pos it ive solutions of trlre 
equat ion 

t-an X n + a X' n = 0 ( 20) 

the first 15 of which are listed for ct. = 1 and a = 5 
in table 1. The expressions for: the spanwise normal 

oh . i ■ ■ " - m 

stress cr_ = and the shear stress t = — -c — follow 

5y 

in the form ' - " ( _ 

°° - - -v - I- \ £_ \ 

i V 1 cog X„ cob \n r: J E a w 

(31) . . 



J 1 \" cos X n coe X n i v E n w 

ja + 1 „-=-!■ i + a c oB 2 X n J 



r-i cos X n sin. X n — «/ E n w 

t ^ s / ~ f- ) ~ r e (22) 

t w 



^='i 1 + a cos 2 X n 
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In the case of an unflanged panel, for which a =, 0 , 
the series "become Fourier series and can "be summed in. 
closed form, giving 



H = 



2F —l / 
^ tan 



tanh— — tan j 



(23) 



T sinh | 



x tw cosh| + cos n 



where 



and 



T= - 



F AT sin n 

tw v E cosh | + . cos T) 



? yi" w 



-n = 



The flange normal and shear stresses, as well as the 
normal stresses along the center line of the" panel, for 
the cases a = 1 and a = 5, are plotted in figure 5(a), 
while the normal stress distributions along ah edge and 
along the center line for an unflanged panel {a =- 0) are 
given in figure 5(h). (See table 2..) It is seen that for 
a = 1 the shear— lag function' cr x (x,w) — a x (x,o) "becomes 
negligible at a distance from the end section equal to 
about 2-| times the width of the panel (x = 5w) , while if 
cc = 5 appreciable shear lag is f present up to a distance 
of about 3 times the width of the panel. If no flanges 
are present, the normal stress becomes infinite at the 
point of load application; at a distance "of about times 
the width of the panel, however, the spanwise normal stress 
becomes practically constant over the cross section. 

Stresses in a Long Antisymmetrically Loaded Flanged Panel 
without Central Stiffeners 
[Panel 2 ] 

If equal and opposite, forces ±J are acting on the 
flanges of a symmetrical panel (fig. 4(b)), the stresses 
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are determined from the stress function 1 

r 

H = li 



I r fz\ 3 _ i 1 

+ i) L W 3 J 



t j 2 (3a 

sin \_fcos Xil — 



{ - '* y sin 



-« ) 



(34) 



where again a = — and the parameters \_ satisfy the 

tw 

equation 

tan X n = ^ (25) 

a \ n 2 + 1 

The first 15 of these parameters 'are listed for a = 0, 
a = 1, and a = 5 in table 3. 

The expression for the spanwise normal and shear, 
stresses in the panel follow from equation (24) by dif-= 
f er en.t iat i on , 



y 

e» sin X_ sin Kf, — 

n n w 



- A 

</ E n W 



a_ = -?-<! Z + 2 ) : e V(26) 



2a+1 W ^. . /Bin \„\ 3 



11=1 1+a sin X n - 



/Bin\ n V 

^ £ sin xjcos X n Z i 'iij^ft) - /f\ *- . 

T. -a/*i) - : ^ ^-^ e ^ W (27) 

^ E t w - — < , . - NS 

1 + a sin 3 \ n — — ^ 

If no flanges are present, the expression for the 
stress .function reduces to 



f — , N 8 - oo COS X„ — i — COB X_ — M: ^-v, 77 "1 

^, 2 -4^ — > 

j 2 |_W 3j x Bln x f 



iThe irrelevant added constants present in equation (24) 
and in the expressions for certain of the other stress 
fiAictions in the following problems are due to the conven- 
tion that H(o,o) « 0. 
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where 

tan X n = \ n """'( 29 ) 

The flange normal stress, as well as the shear ; stress 
at the flange and along the center line, is plotted for 
a = 1 and a = 5 in figure 6(a), while for an uhfTanged 
panel the normal stress along the edge and the shear~~s~tr es s 
at the center line are shown in figure 6(h). (See tahl e 4. 
It is seen that in all cases the normal stresses closely 
approach their limiting values at a distance from the end 
which is less than the width of the panel. 



Stresses in a Long Hanged Panel without Central Stiff ener 
Suhjected to Asymmetrical Concentrated "Loadings 



If unequal loads are acting on the flanges of a sym- 
metrical panel, the corresponding stress function can he 
ohtained hy superimposing proper multiples^ of the symmet- 
rical and ant isymmetr ical stress funct ions - (l8 ) and (24). 

For example, if, as in panel 3, a force F is act- 
ing upon the. flange along the edge y = w, the other 
flange "being unloaded (fig. 4(c)), the relevant stress 
function is given hy one— half the sum of the functions 
(18) and (24). The corresponding stress distr ihut ion is 
indicated for a = 1 in figure 7. 



Stesses in a Long Symmetrical- Flanged Panel with 

Central Stiffener 

C Panels 4 and 5] 

If arhitrary concentrated, loads are introduced into 
the two flanges and into the central stiffener of a long 
symmetrical panel, the loading may he considered as Tihe 
superposition of two loadings:" (l) equal loads F 6 act- 
ing on the two flanges and a load F^ acting "on the stiff- 
ener, and (2) equal and opposite loads ±F introduced into 
the flanges, no load acting on: the stiffener. In the sec- 
ond case (fig. 4(e)), the stiffener is ineffective and_the 
stress distr ihut ion is ohtained from equations (247 to~ 
(27)„. That the stiffenex is Ineffective follows from The' 
convention that the stiffener is concentrated along a 
line. If the width. of the stiffener is small in compari— 
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son with, the width of the panel, the effects actually 
present would in any- case he negligible. 

In the case of a symmetrical loading (fig, 4(d)) con- 
sideration of the shearing action, along the two edges, of 
the central striff-ener indicates that the shear stress is 
discontinuous at the stiff ener. As a result, it is found 
that the stress function has different analytical expres- 
si-ons-in the two sections of— the "panel . Because of- the 
symmetry, however, it is suf f icie:nt to determine the state 
of stress in one— half the panel. ' In the positive half of 
the panel (o < y < w) the stress, function is obtained in 
the form 



H = 



t(a + a m + 1) 



OD 

_ T 



n=i 



r 



2 



*r i 
n ^ 



sin X n - + a m 



cos Xjj- — « 



— X — 
E n w 



(30) 



where the. constants are defined as 



1 n 



w 



A 



a = 



e 
A 



m 



2tw 



J 



sin \ n — + a m cos 



dy 



(31) 



(32) 



K n 88 




a 



+ 3" 



+ am )} { ^ 



s in \, 



cos \ 



•n 



£1 

wj 



dy (33) 



and the parameters 
equation 



a \. n tan \ n — 



are the positive solutions of the 
tan \„ 



a + On 



(34) 
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The stresses a x and t follow from equation (30) 
"by differentiation in accordance with equation (7). 
Although the expression in equation (30) is somewhat com- 
plicated, it is found that, except in the immediate .neigh- 
borhood of the end x = 0 where the stresses are already 
known, the presence of the exponential factor in the series 
"brings about such rapid convergence in practical cases 
that very few terms of the. series are needed for -ordinary 
accuracy. In a later section the stress distr ibut ion 
associated with equation (30) is evaluated numerically for 
the case of an actual test panel (panel 10) considered pre- 
viously in reference 3. ' ~ 

Panels 6 to 9 represent a special case in wnich 
Aja = 2 A e . If the area of the central stiffened is just 
twice the area of' each flange, the expressions involved in 
the solution are considerably simplified." Tn this case a 
general symmetrical loading can he treated^ by superimpos- 
ing two basic loadings: (l) eq*ual loads 3P acting on the 
flanges and a load 23f acting on the stiff ener (panel 6), 
and (2) equal loads F acting on the flanges' and a. load 
— 23P acting on the stiffener (panel 7).. Again, because erf 
symmetry, it is sufficient in each, case" to consider" the 
positive half of the panel. 

Tot panel 6 (fig. 4(f)), each half of the panel behaves 
in the same way as panel 1. She stress function is obtained 

by replacing w by and y by ^2 ^ — 1^ in the ex- 

pression for the stress function for the complete paneli 
so that, for o < y < w, * 

^ cos L sin X (.^r, — l) 



t|2a+lV w J n^.* n(l .j +3a C os a \ n ) 

ff = JU-L^ 4 - — — V " 7 e V* " ■ w ' f- (36) 

" tw | 2a+1 1 + 2a cosfe X n J 



H 35) 



and 



/ , r V ' /i 



T = 4 /£ £- >. cos K * sin bi! I ~ V e" 2 ^ 1 ^ * (37) 



x 



B tw 



n=i 1 *" 2~a "cos* K n 



A 



e 



where a = and the constants X are now the positive 

t w 

solutions of the equation 
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tan X n + 2a >v n = 0 (38) 

The stress function corresponding to panel 7 (fig. 

F 

4(g)) is obtained by using 1 e = - ~B- = J? and A m = 2 A e 

in equation (30). After some simplification there follows 
for o < y < w, 

t ^ x n(l + 2a sin 2 X n ) J 

"2, / y ^ x 

■p \ s in X n s in \ n ( 2 - 1 : — 2 /= \ n — 

*_ = 4 — L ~ ^^-J* 1 0 v . n w (40) 

tw n=i L sin 2 x. 



n 



and °1 / y \ 

_ V sin x n cos x n I 2 - - 1 ! 

/g jr_ > n n v w y 

■n/E tW n=1 - - , 2 . 



' -2 /- X n - 



t-4./5^: £ ^— VE (41) 

1 + 2a sin® X a 
where a 



t w 



and 



cot X n - 2a \ n = 0 (42) 

The first 15 of the solutions of equation (42), for 

1 5 
a = — and a = are given in table 5. 

The spanwise normal stresses and t"he shear stresses 
at a flange and at a quarter section (^y = -gy are plotted 

for panel 6 in figure 8(a) and in tablo 6, ■ for a = 
and a = Jy, while in figure 8(b) and in table 7 are pre- 
sented curves representing the flange normal and shear 
stresses for panel 7. It"~ris to be iremarked that- ths 
stress cr x is symmetrical about the center line (y = 0) 
in both casos, symmetrical about the quarter section 
f w 

y = ~o \ for panel 6 and ant isymmetr ical about that line 

V J 

for panel 7. 

By combining suitable multiples of the stress func- 
tions given in equations (24), (36-), and (39), solutions 
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may be obtained for: an arbitrary set of three concentrated 
loads introduced into the flanges and the stiffener for 
the case in which A m - = 2 A e . In particular, if equal 
loads 3? are acting on the flanges, while no load acts 
on the stiffener (panel 8, fig. 4(h;), the r elevantf" slir ess 
function is given by one— half the sum of the functions 
given in equations (36) and (39). The normal stresses in 
a flange, in the stiffener, and at a quarter section, and 
the flange and stiffener shear stresses are presented, 

for a. = -i and a. = ~ t in figures 9(a) and 9(b), re spec— 
2 <2 

tively, and in table .8. If a load 2P is acting on the 
stiffener, while no loads are introduced into the flanges 
(fig. 4( i ))» the stress function is given by one— half the 
difference between the functions of equations (36) and 
(39). In the present case, however, it is seen that the 
expressions for the flange and stiffener stresses, respec- 
tively, in panel 9, are obtained by interchanging the cor- 
responding expressions for panel 8. 

For the purpose of comparison there are included in 
figure 10 curves representing the normal stresses along 
the longitudinal edges of unflanged panels subjected to 
(l) equal concentrated axial forces 3T acting at the cor- 
ners of the free edge and a force 23? acting at the cen- 
ter of that edge, (2) equal forces P ' acting at the cor- 
ners and a force —23? at the center, and (*3) equal and 
opposite forces ±3? acting at the corners. (See table 9 

Panels 6 and 7 can be considered sections of a semi- 
infinite sheet stiffened by a series of equally spaced 
longitudinals of area 2A e and constant separation w. 
where, for panel 6, equal axial loads 23? are acting on 
successive stiffeners, while, for panel 7, successive 
loadings are of magnitude 23? but are alternately com- 
pressive and tensile. By superposition the stress dis— 
tribution is obtained for an arbitrary pair of cohceh— ~ 
trated loads repeated periodically along the edge of the 
sheet. 3Tor example, the stresses given for panel 8 ( or 
panel 9) can be interpreted as the stresses in a section" 
of a sheet of this type where alternate stiffeners are 
not loaded. Similar statements apply to the curves pre- 
sented in figure 10. 
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Stresses in a Doubly Symme'tr ical. K-ect angular Bos Beam 

Clamped at~One End and Subjected to Concentrated 

Bending loads Applied at .the Free End 

The rectangular box beam shown in figure 2(a) is sub- 
jected to concentrated "bending loads applied at the free 
end. It is here supposed that no longitudinal, center 
st iff oner is present, although the analysis of a later 
section includes also the theory in the more general case, 

The stress function is determined in the form 



r 



H = wa 0 ' 



where c 
stress , 



_ f^cos \ n sin \ n - sinh |j, n - 

f 2W X )/5» l-r- 1 — 



C-lnSlK n c (i+ ^ cos X n ) coah p 
is the maximum elementary spanwise normal 



k (43) 



M( I) h _ Fjh 
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'(44) 



P is a dim ens ionles s paramet-er representing the ratio 
of the moments of inertia of the web and the sheet about 
the transverse principal axis of the beam and is given as 



(45) 



the parameters 
equation ^ 



n 



are the positive solutions of the 



tan X n + & X n = -0 



n 



(46) 



which is of the same form as equation (20), and is 
defined by the relationship 



M-n 



^ 1 x 
E w n 



-(47) 



The stresses a x and T follow from equation (43) 
by dif f-er ent iation 
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CT x = CT o 




CO 



f + acp+i) Jl ~ ) — — ^ — 

l V & l i_. \ n (l + P cos 2 \ n ) cosh n n 

11=1 



cos \ n cos \ n I sinh p, n 



(48) 



t = -a, 



1 w n=i x n (l+P 



^f" cos X- sin X„ — cosh 



w 



^n 



cos 2 Xn.) cosh (j, n 



(49) 



It is seen that the series parts of equat i ons ( 48 ) and 
(49) are of the nature of corrections to the linear 
stresses predicted by the elementary theory. lor span- 
width ratios in the practical range the parameters |j, n " 
are of such a magnitude that, except in the immediate 
neighborhood of the fixed end (x = l), sufficient accu- 
racy is obtained "by retaining only a small number of the 
terms in the series. 

An alternative expansion for the stress f unct io n _can 
he obtained by conventional methods, which are somewhat 
simpler than those employed here," in the form 



CO 



(-1) 



n+i 



H — 5" CT 0 ' 

tt I s — (2rL—l)~ sinh8 n + B 8 n cash 8 n 



sinh 6 >, — , . 

sin U-zlk|(50) 

2v I. 



where 



8^ = 



(2n - 1) 



2 
I. 



TT AS W 



(51) 



■p = 



from which there is obtained 



CO 

V 



4 I £e w _ 
n TOG T 0 " 



/ ., \ n+.i cosh 6-n — , _ ■ " v ' 

; ^ — ^ i 2n - l)TT 5 (52) 

n=i 2h— 1 sinhS n + BS n coshS n 2 I 



u% t **<=i£i 6lah 6a .» ... i^ik 5 (63) 

rr I s I 2n — 1 sinh 8 a +£ 8^ dosh 8 n 2 I 



This form of the solution has the disadvantage that the 
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series, instead of being of the nature of end corrections, 
serve to determine completely the stresses. Moreover, 
the oonvergonae is least rapid in the neighborhood of--the 
flanges (y = ±w) where the magnitude of the stresses is 
generally of greatest interest. Finally, the procedure 
leading to equation (50) cannot be directly employed in 
dealing with beams -with cut— outs while, as will be shown, 
the present methods apply without modification. 

The normal stresses in the flanges and along the 
center line of the sheet for a beam with length— width 

ratio — = 2,5 are shown as the solid curves in figure 
2w 

11(a) for 3 = 1 and in figure 11(b) for the limiting 
case 3 = 0. (See table 10.) The solutions obtained by 
Eeissner and by Kuhn and Chiarito for the former case will 
later be compared with the present results. 

The transverse distribution of normal stress at the 
root- is presented for 3-1 and 6 = 0 in figure 12 
and in table 11, The parabolic distribution predicted by 
Eeissner' s solution for the case 3=1 is shown as a 
broken— line curve. The fact that : the stresses in the 
flanges and in the central fibers of— the beam, as pre- 
dicted by the Eeissner solution and by the present solu- 
tion, agree almost exactly along the entir-e span except . 
for small deviations at the root, as is shown in a later 
section, indicates that the parabolic approximation is 
extremely accurate except in the immediate vicinity of the 
root. In the limiting case 3=0 theoretically infinite 
stresses are predicted at the root corners. 

The so— called effective width of the cover sheets, 
definod by the relationship 

w 



can bo expressed as ■ fallows in terms of the flange stress 
and the elementary linear stress 




(54) 



x 



w 



off 



w 



= 1 - ( 3 + 1) 1 



v x { x,w; 



(55) 



The solid curves of figure 13 represent — the ratio of the 
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effective width to the actual width 
a function of distance from the fre 
P = 0. (See tattle 12;) For p = 1 
Width is found to be 78.6 percent o 
compared . with the value of 80.1 per 
Eeissner. In the limiting case p 
90 percent efficient at distances' f 
than one— fourth the width of the sh 
width rapidly decreases to zero as 



of the cover sheet as 
e end' for P = 1 and 
, the effective root 
f the actual width as 
cent obtained by 
= 0, the sheet is over 
rom the root greater 
est, but the effective 
the root is approached. 



Effect of a Cut— Out on the Stress Distribution 

in a E-ectangular Box Beam 

Subjected, to Concentrated Free— End Bending Loads 

If a section of the cover sheet is removed at x = I 
(fig. 2(b)) and no sheet stiffener is present, the stress 
distribution in that part of the sheet between the free 
end and the cut— out is derived from the stress function 



co y x 

V* cos X n sin \n ~ s inh u,„ t 
. x y > w 11 I 

K= v cr 0 ^ T - + 2(p+l) L — : 

1 v X n (l+P cos \ n ) sinh ti n 

where the constants are defined as in equation (43). 
follows that 



(56) 



°.)7 



+ 2(P+1) 



^ cos X n cos \ n - sinh n n T 



It 



(57) 



n= 1 1 + p co 6 ' 



s inh 



t= — o , 



1 w v E 



+ 8(p + l)/Ji r 

w v E w i_+ p cos 2 



^— cos X n sin — cosh jj, n — 



sinh n n 



> (58) 



The spanwise normal stress distributions in a flange 
and along the central line of tihe sheet for the span— width 
I 

ratio — = 2.5 are shown as broken— line curves in fig— 

2w — 

ure 11 for P - 1 and P = 0, while the variation of the 
effective width along the span as computed from equation 
(55) is presented for P = 1 and p = 0 in figure 13. 
(See tables 12 and 13.) A discussion of the results is 
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postponed until a later section in which a comparison is 
made with the solution given by the procedure of reference 2, 

Stress Distribution in a Semi-Infinite Sheet with One Stiffener 

If a load P is introduced through a stiffener of _ 
cross section A into the straight— edge of an infinite 
half— sheet assumed rigid in the direction parallel to the 
straight—edge (fig. 3(a)), the. stress distribution is deriv- 
able from the stress function 



where 




.(59) 



(60) 



so that, in particular, the expressions for the normal 
stress in the direction of the stiffener and for the shear 
stress take the form 



21" n ~7 P sin k p x , 
o- x= f e — (61) 

P + -X" 



2kJ f° ~ y 5 cos k p x , . 

T- — / e d p (62) 

0 P + T" 

Along the stiffener (y = 0) these expressions can be writ~ 
ten in terras of tabulated functions as 

, \ 21 f 2ktxfn „. 2ktx 1 2ktx 2ktx\ , . 

. ■ 2kJ r 2ktxfTT _. 2ktxT 2ktx rt . 2ktx\ . 

t(x,o) = — -ssin — - — - - Si — - — - cos* — - — Ci — - — > (64) 
■ * nA » v A {2 A J A A J 



with the conventional abbreviations 
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Si' x = / 



sin u du 
u 



Ci x = 



cos u du 
u 



(65) 



J 



In the limiting case A = 0, when the sheet is un— 
stiffened", » 

j f — yp sin k x P ? — l kx 



H = — 



TTt 



e 



P 



dp = tan 

TTt 



(66) 



F ! kx 
CTx " ^t y 2 + k 2 x 2 



kg y 

nt y z + k 3 x 2 



(67) 



and, in particular, the normal stress in the direction of 
the applied load is given "by . 



a x (x,o) 



TTtkX 



(68) 



The corresponding stresses in an unstiffehed isotropic^ 
aheot; are known to. "be of the form 



2Ex 



T = 



TTt(x 2 + yB)" 

2F.x 2 y 
TTt(x 3 + y 2 ) 2 

22 1 



TTtX 



(69) 



It seems of some interest to compare the stress dis- 
tributions just obtained with that associated with a load 
uniformly distributed over a snlall area at the margin of 
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tho shoot, " In particular, if the load F is uniformly 
distributed over an area t 2 (fig, 3(b)), tho stress 
function takes the form : - 



co 



P f -sxb sin e* .. ■ 
H = ^ra~/ o _ sin py dp (70) 

0 

and tho normal stross along the lino of symmotry follows 
by differentiation, 

, V OK (X,o ) P; •, "I / t \ . f„. V 

c,,(x,o) = — 2 — = ~s tan ( — ) (71) 
oy tt t V-^/ 

The corresponding stress for an isotropic material is 
givon by the expression 

x 

*«<*.«> - ^ * * r h^t t72) 



In figure!4 the stiffenor shear and normal strossos 
are represented in their depondonce upon trho ratios 

FGc A F 1 ' 

— — and — , whilo in figure 15 the flange normal . 

t 3 A 

stresses corresponding to equati-ons (71),' (72), and (63) 
are compared. (Soo tables 14 and 15.) 



MATHEMATICAL FORMULAE 101? OP THS B OXJISTD AR Y— V ALU3D PROBLEMS 
Strossos in Axially loaded Symmetrical Panols 

Symmotr i.eal loading.— The notations and dimensions 
assumed in this section are indicated in figuro l(a). 

Bocausc of tho symmetry tho transverse displacement 
aust vanish along tho lino y = 0 and hence, from equa- 
tion (4), must vanish everywhere. : Thus tho differential 
oqu'ation (8) for tho stross function E tafcos tho form 

^H-^i^='o (73) 
ox 2 E dy s 



At tho froo ond ( x = 0) the spanwiso. normal stross 
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a-j. = must vanish excopt at ^ = 0 and at %J = ±w 

and Is finite at those points, so that H(o,y) must "bo 
constant. Sinco, according to eqjiiation (?), the stresses 
remain unchanged if H is increased by a constant, it is 
permissible to require that 

H(o,y) = 0 (74) 

Because of continuity of tho spanwiso strain e x , 
it follows that tho spanwiso stress a x = ^ is continu- 
ous at y = 0 and tho condition of equilibrium between 
the sheet and stiffenor along the lino y = 0 can be 
written in the form ... 

A H — ^r— ! = -t [t(x,o+0 -t(x,o-)] (75) 

ox 

If equation (75) is written in terms of the stress func- 
tion H and if the fact that T is ah odd function of 
y is used, this condition becomes 



OS 



HCx,o4 - &H(z ' o) 1= 0 
2t . oy J 



r (76) 



s A m dH(x,o) 

Hv.s,o+) — — = constant 

2t by f 

It now follows from the end condition 

A«-a x (o..o) ^ i« m H ,0) = ? m (77) 

and from equation (74) that tho constant in equation (76) 

F 

m 

has tho value of — ~, and the condition along the pos— 

ft V 

itive side of the line y = o takes the form 



(78) 
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Since -H and -l^- aro odd functions of- y, it is soon 
■ Ox 

that thoy aro discontinuous at the lino y = 0 (except 
at x = O) unless A m = T m = 0.' Also, from oquations 
(74). and (78), it follows that 

~> 

oH( 0> y) „ n 



1 im 

dH( x,o) J m 



lim 



by 



A 



m 



(79) 



The xjonditT.on along the l.ine y = w can bo obtained 
from the spanwise condition of equilibrium 

v; 



St 



J cr x dy + 2 A e a^x.w) + A m .a x (x,o) = 2:F 0 +]r m (80) 



If a >: - ; is introduced into equation (80) and if 



oquation (78) is taken into account, this condition 
roducos to 



dH(x,w) 
t H(x,w) + A c — ^ = J { 



(81) 



At the end x = I a condition of the form 



a 0 H(l,y) + a x 5H ^' y) " ? f(y) 



y > 0 



(82) 



is assumed. Among the s itmations ■ in which equation (82) 
applies, tho following cases may be notod: 

1. If the end x = I is clamped so that the span- 
wise displacement u vanishes, it follows from equation 
(5b) that 



SH( l,y) 
Ox 



( 82a) 



2. If a section of tho panel, including the central 
stiff enor, is removed at x = I, it follows that 
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cr v( 1 »y) = O if y ^ iw so that 'H(l.y) = constant. 



Equation. (78) then requires that H(l,o) 
condition at x = I "becomes 



S and the 

3t 



y > o 



It should be not iced that , since 
of y 

+ 3F 

HU.y) = 



m 



3 t 



(82b) 

is an odd function 



y < 0 



so that H(l,y) is discontinuous at y = 0 unless 3? m = 0 

If F m =0, it follows from symmetry that the panel can 

be considered as clamped along the line x = ^, so that 

. 2 - 

this case is then reduced to case (1). 

3. If a section of the sheet is removed at x = I 
but tne central stiffen e r is "built in, it again follows. 

H(l,y) = c 



that 



y > 0 



1 



(82c) 



Although constant c cannot in this case be determined 
from static considerations, if- the stress function is 
determined in terms of c, the additional condition 

u(l ,w) = u.(l ,o) 

which becomes, from equation (56), 

w 



ox 



dy = 0 



is sufficient for the determination of the constant. 
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4, If the end x = I is Subjected t o an axial load- 
ing distribution t g(y) vrhich. is statically equivalent 
to the loading at" the end x = 6 -, eo that - 



a (I ,y) = — = 

&y 

it follows from equation (78) that 



g(y) 



y 



r 

H(l ,y) = - j g(u)du - 

o 

y > 0 



•"-m 
2t 



(o) - -a 

2t 



(82d) 



If the- preceding results are summarized, the stress 
function _H is completely determined in the region 
[o < x < J, o<y<w] and hence, by symmetry, over the 
pntire panel, by the following s-et of —e~q.ua t ions : 



ox s E dy 3 



H(o,y) = 0 

a 0 H(I,y) + ai «iil?L> = f (y.i 

ox 

2tH(x,o + ) ,A m ^|^ = -F ffi 
tH(x.w) + . A e = F e 

5y 



<A1) 
(A2) 
(A3) 
(A4) 
(A5) 



Ant i s.vmmet r ical loading,— The notations and dimen- 
sions assumed in this section art indicatfid in figure 1(b). 

Since the transverse displacement v does not van- 
ish in this case, the stress function is governed by the 
more general dif f er enfrial equation (8) which can be writ- 
ten in the form 



by\ dx s B dy 3 / ~ 



(83) 
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The conditions along the lines x = o, x = I, and 
y = w are identical with those in the preceding section" 
and are given in equations (74), (81), and (82). Since 
ct x (x,o) = 0 from symmetry, equation (78) is replaced "by 
the condition 

SSktSLl . 0 ' (84) 
oy 

An additional condition .is obtained from the equa- 
tion of moment equilibrium which can be written in the 
form 

r 

t / y ct x dy + w A cr x (x,w) - w P (85) 

0 

cH , . 

If a_ = — — is introduced into equation (85) and the first 
x oy 

term is integrated by parts , - 

w 



t H(i,r) - J Je dy + A = » 



(86) 



and equations (8l) and (86) are compatible only if H 
satisfies the additional restriction 



J 



H(x,y) dy = 0 ... (87) 



The following set of equations then determines the 

stress function H in the region £0 < x < 1 » °. < 7" < w 3 : 

_&./&!5 + a = 0 CBi) 

dy Vox 2 E oy 3 / 

H(o,y) = 0 (B2) 

a 0 H(l,y) + ai SH( *' y) = f (y) ( B3 ) 

d H (* ? °I =0 (B4) 
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„ / \ dH(x,w) 
t H(x,w) + A = I 

£>y 



(B5) 



j H(x,y)dy = 0 



(B6) 



Stresses in the Cover Sheets of Doubly Symmetrical Box Beams 

The notations and dimensions assumed in this section 
are indicated in figure 2. 

If it is assumed that the "bending loads are applied 
symmetrically about the line y = 0, it follows as before 
that the transverse displacement ;v vanishes identically 
and the stress function is determined by the differential 
equat ion (73 ) . 

The conditions along the lings x = Q and x = I 
are identical with t"hose derived in the preceding sections: 
namely , 



H(o,y) = 0 

a 0 HU.y) + ai 2Elhll = f (y) 

ox 



(74) 
(83) 



where equation (82) includes, in particular, the case when 
the beam is completely clamped at x = l' [a 0 = 0, f (y) = O] 
and the case when the side structure is fixed but the 
sheet and 1 ongitudinal are not attached to the support 
[a x = 0, f (y) = 0] . 

Along the positive side of the line y =» 0 , the dis- 
cussion of the preceding section again leads to equation 
(76). In the present case, however, the normal stress 

dH - 
cr y = vanishes at x = 0 for all y, and it fnorllows 

from equation ('74) that the c ons t ant~ of integration in 
equation (76) must vanish, giving the condition 



(68) 
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The equation of moment equilibrium about the trans- 
verse principal beam axis can be written 



w 

I 



2th f cr x dy + — ^ a x (x,w) + h A m a x (x,o) = M(x) (89) 

} ■ h/2 
o 

where M(x) is the moment of the applied- bending loads. 
If equation (89) is expressed in terms of H and if 
equation (88) is taken into account, the boundary condi- 
tion along the line y - v follows in the form 

£b h(X|W ) + J ' M^l" = n'(x) § (90) 

w oy 2 

In summarizing these results, it may be seen that 

the stress function E is determined in the region 

[o < x < I, o < y < w] by the ' f ollowing set of eq.ua— 

t i ons : ■ ... - — — 

&!h + ff £3 = o ' (ci) 

• ox 2 ' E dy* ■ 

H(o,y) = 0 (03) 

a 0 H(l ,y) + a x i=iL*li = f (y) (08) 

ox - 

2t H(x,o+)- A m SHC:: ' 0) = 0 (04) 

oy 

H(x,w) + I„ = M(x) I (05) 

w oy 2 

SOLUTION OP THE PEOBLEMS 



One of the difficulties encountered in solving the 
sets of equations derived in the preceding section arises 
from the nonhomogene ity of the conditions at both the 
boundaries x - constant and ' y = constant. In most cases 
however, the problems considered can be reduced to prob- 
lems of a more conventional type by the definition . .. . 

H = $ + P (91) 



where P is a particular solution of the governing differ 
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ontial equation (73) . or (83), determined so as to satisfy- 
either the end conditions along the "boundaries % s=. con- 
stant or_ the oquilibr ium conditions prescribed along the 
boundaries y = constant, so that the corresponding con- 
ditions become homogeneous in $ .. In case "both procedures 
are possible it is pr ef"5rr able to. satisfy tho equi lihrium 
conditions for tTtrreasons: (l) The function P so 
determined is then a stress function corresponding to the 
state of stress in a sheet subjected to external forces 
that are statically equivalent to the actual forces and, 
according to Saint Tenant's principle, if the length— width 
ratio of the sheet is sufficiently large, tho two stress 
distributions will he nearly identical except in the 
regions near the ends of— the sheet- Thus the function $ , 
which must "be expressed, in general, as an infinite series, 
is of the nature of an end correction. If, however, the 
particular solution P satisfies only the end conditions 
along the. . s ect i ons X = cons t ant— the infinite series 
representing $ will bo in general significant over the 
entire span. (2) Certain computational advantages follow 
from the fact "that if- -P is made to satisfy the equilibr-ium 
conditions the variation of- - $ ih the spanwise direction 
is expressed in terms of exponential (or hyperbolic) func- 
tions, while in the alternate procedure the reverse is 
true. If the condition of support at x = I is such that 
f(y) - 0 in equation (82), as is true in many practical 
cases, the first procedure is applicable with P = 0. Cer- 
tain special cases of— the problems considered here have 
"been solved by this last method (references 3 and 4), hut 
the solutions obtained in this way ar e not so well adapted 
to numorical computation as thoae obtained here. 

In the rest of this section the boundary value prob- 
lems formulated in the preceding section are reduced to 
sets of ordinary differential equations with terminal con- 
ditions, all of which can be solved "by elementary methods * 



The Symmetrically Loaded Panel 

A solution of equation (A) is assumed in the-form 

H = Z X n U) Y n (y) +-p(y) (91a) 
n 

Now the differential equation (Al) requires that X n , Y n 
and p satisfy the ordinary differential equations 



P"(y) = o 



(92) 
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(93) 
(94) 



where X n and are arbitrary constants s at isf y ing . t he 

relet ion 



(95) 



The boundary conditions of equations (A2 ) to (A5) take the 
form 



n 



J 



*n (y). 


= -p(y) 


(9.6 a) 


T n (7) 


= f (y) - a 0 p (y) 


(96b) 


X n (x) 


= -j2t p(o) - A m p 


(o) + fJ (96c) 




= -|t -p(w) +■ A e p 1 


(w)-» e l (96dL) 



The equilibrium conditions epilations (96c) and (96d) 
are satisfied and homogeneous conditions- of" the type 
required' are obtained if it is required that 



and 



2t P (o) - A m p»(o) = - P m 
t p (w) + A e p' (w) = P e J 

2t T n (o) - A m Y n ' (o) = 0 
t Y n (w) + A e Y n '(w) = 0 



(97) 



(98) 



Equations (92) and (97) are sufficient to determine the 
particular solution . p (y) , while the differential equa- 
tion (94) and the homogeneous boundary conditions (98) 
constitute a "characteristic value" problem of the Sturm— 
Liouville type over the interval 0< y< w. The constants 



'n 



are determined as a discrete infinite set of charac- 



teristic numbers; and, corresponding to each such \ 



34 



HACA "Technical Note Vo. 894 



a characteristic function Y a satisfying the conditions 
of the problem is determined. It is a property of such 
systems that the set of characteristic functions so ob- 
tained is orthogonal over the interval involved and that 
any function having a continuous second derivative in 
that interval has a unique expansion in a series of those 
functions. ffur thermore., if the. permissible oondition of 
normalization 



T n s dy = 1 (99) 



o 

is imposed, the constants in the expansion 

T(y) = S c n Y n (y)~ 
0 < y < w 
are determined by the formula 

■xV 



(100) 



I-(u) Y a (u) du (101) 
o 

With p(y) determined and tho relationships neces- 
sary for the determination of 'the orfrhonormal set 

• Y n (y)>- known, it remains to determine X n (x) so that 
I . J 

equations (96a) and (96b) are satisfied. If it is assumed 
that f (y)- ha's a continuous second derivative in (o,w)', 
it follows from equations (96a) and (96b) that, according 
to equation (101) X n (x) must satisfy the conditions 

w 

X n (0) = - T p(u) Y n (u) du (102) 
•-'o 

w 

a D X n (i) + * x X n '(i) = J-^tM - a 0 p(u)| Y Q (u)du (103 ) 

and equations (102) and (103), together with equations 
(95) and (93), are sufficient to determine X n (x). 

The solution of the boundary value problem of equa- 
tion (A) in the region [b < x < I , o < 7 < v] is given, 
as a summary, by 

H = g X n (x) Y n (y) + p(y) 

where 
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p" = 0 

St p(o) - A m p' (o) = -T m 
t p (w) + A e p r (w) = ! e j 



Y a = 0 



2t Y a (o) — A m Y n ' (o) 



0 

t Y n (w) + A e Y n *(w) = 0 



w 



dy = 1 



X n (o) 



1 



a Q X n O) + ai X n '(l) = f[ (f - a 0 p)Y n dy 

o 

Explicit solutions in several cases are presented in 
a preceding section. The determination of H is illus- 
trated "by considering the first case in that section ' 
(panel 1 ) . 



Here 



1 = CO 



*e = * 



= f (y) = A m = P m = 0 



A e = A 



J 



(104) 



The particular solution is ohtained as 

P (y) = cj y + c 3 
p(o) = 0 t p(w) + A p * (w) =3* 



P(y) = 



t (a + 1) 



a = 



tw 



(105) 



£6 
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Next the functions Y n (y) satisfy the conditions 
Y n (y) = ..Oj sin \ n £ + c 2 cos \ n - 



t r n ( w ) + A Y n '(w) = 0 



w 



dy = 1 



from which the characteristic numbers \ n are determined 
as the. positive solutions of the equation 

tan X Q + a X n = 0 (106) 

with corresponding normalized characteristic functions 

y 



Y n (y) = y n sin X n J 



w(l + a cos 2 k n ) 



(107) 



The terminal, conditions governing the function X n (x), 
which according to equation (93) is of the form 



v 3 n w . „ „ VE 



c i e 
(102) and (103.) 



~ X n * 

+ c 3 e v E w , f oil o w from equations 



X n »(«0 = 0 



w 



, v ' n A p „ v wY n 
X n (o) = — £ sin X n i dy = 

n t-(o + i) 7 w w t 



W Y n J 1 COS X, 



(108) 



from- which ,tr 

._ /G x x 

v ^ N _ w Y n I* -cos X n «/ E n w 

A n (, x ; - - — ■ — — e 

t \ n 

and the solution is obtained in the form 



1 tf r 

"+ 2 / 



f 2 ! cos h n s in \ n £ 



w 



8 *n 2 
3 n W 



a 4 ^ X n (l+a cos%) 



(109) 



(110) 
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where 

A. 

tw 



and 



tan K n + a \ n = 0 



The Tcliiity of this solution can "be rigorously 
established if the fact is taken into account that, for 
large values of n , 

x m IS* - Hz 

n 2 

(-) n+1 2 
cos- \ n « — 



a(2n - l)rr 

The Antisymmetrically Loaded Panel 

A solution of equation (B ) is assumed in the form 

K = 2 X n (x) Y n (y) + p(y) - -(Hi) 

n 

where, in consequence of equation (Bl) , the functions 
X n , Y n , and p must satisfy the differential equations 



P"* (y) =o (112) 

X/U) -(t 1 ) 3 X n (x) = 0 (113) 

V"<y) + (if) 3 Y n'(y) = ° . . < 114 5 



and, as "before, 

Equations (32) to (B6) are now satisfied if it is 
required that 

gX n (o) T n (y) = - p(y) (116a) 

\ \ a o X n< l > + a i s n'( l > r Y n(y) = f (y) - a 0 p(y) (H6h) 



38 



IS ACL Tochnical Note- No. 894 



p» Cp) 



t p (w) 4 A p ' (w) = F 



p dy = O 



(117) 



Y n '(o) = 0 
t Y n (w) + A Y n «(w) = 0 



Y ft dy = c 



(118) 



The conditions given ,in equation (117) together with 
the differential equation (112) uniquely determine p(y), 
while Y n (y) must satisfy the boundary value problem con- 
sisting of the differential equation (108) and the homo- 
geneous boundary conditions of equation (118). Although 
this problem is not of the conventional Sturm— Liouvillo 
type, it can be shown that the characteristic functions 



Y n are again _ or thogonal 



function F (y ) 
satisfying the 



having a 
c ondit i on 



over _(o,w) and that an arbitrary 
continuous second derivative and 
/ w 3'(y)dy = 0 has a unique 



expans l on 
f unct ions 



of 
Y 



n 



the type of equation 
are normalized, 



(100) where , if- the 



/ 



dy 



(119) 



the constants are given by the formula of (101). Thus, if 
it is assumed that f (y) satisfies the conditions men- 
tioned, equations (116) imply equations (102) and (103). 

The solution of the boundary value problem of equation 
(B) in the entire region [ o < x .< I , [ y | < w] is given , 
as a summary, by 



where 



H = £ X n (x) Y n (y) + p_(y) 

pi" = 0 
p'(o) = 0 
t p(w) + A p ' (w) = 1 

J p dy = 0 
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T n «(o) = 0 
* Y n (w) + A T n i (w) = 0 



vr 

X 



Y dy = 0 



' V -|(^Yx n ' = 0 

•3 V w J 



ay 



x n ( 0 ) = - f P r n 

w 

a 0 X n (l) + &1 X n '(0 = | (f - a 0 p)Y t 



The explicit solution for the case J = oo , a Q = f (y) = 0 
is given in equation (34). 

The Rectangular Box Beam 
If the solution of equation (0) is assumed in the form 
H = Z X n (x) Y n (jr) + P(z,y) (120) 



where 



and 



ox 3 E dy 2 



I 3 X » - p.,, 3 X, 



w s Y» + Y 



2 

n -n 



0 
0 



(121) 

(122a) 
(122b) 

(122c) 
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equations (02) to (C5) take 'the form 

i 

Z X n (o)Y n (y) = - P(o,y) (123a) 
S|a 0 X n (l) + ai X n 'CO} Y n Cy)'=f (y)-ta 0 P ( I ,y ) + a l P 3C ( I ,7 ) ] _.(133T) ) 
S [2t Y n (o) - A m Y n '(o)| X n (x) = -[2t P(x,o)-A m P y (x.o)] (123c) 

= M(x) § ~ [~ +' I w V x ' w) ] < 123d ) 

It is first shown that, if M(x) can he expressed 
as a polynomial in x, a function P(x,y) satisfying 
equation (121) can he determined as a polynomial in x 
and y so that the right-hand sides of equations (123c) 
and (123d) vanish, and 

P( x ,o) - P_(x',o) = 0 (124a) 
2t ^ 

— P(x,w) + I w P v (x,w) = K(x) I (124b) 
It can be d.irectly vor.ified that the expression 

„(,.„ . *{.« [, + ^ § «■«[(£ ♦ £ £> -(-n)J 

* (f >" « l( 5 * £ 5x). - (S + £ 5) * (* * S)] 

- (iY m ti (x) ,7 c + ^ z^-k, * is z±V ks Cz!+f5 zTs 

\aj L\ 7 ! 2t 6.'y \5.' 2* y \3J 2t 2J/ 

-*.(' «-&-)]♦' •} (12S) 



where the law of format ion of the following terms is obvious, 
satisfies the differential equation- (121) and the condi- 
tion of equation (124a) for an arbitrary choice of the 
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constants 



■l > 



-2 » 



The factor 



h 

21 



where 



is the total moment of inertia of the cross section about 
the transverse principal axis" and is given 'as 



i = i + i + h 2 

s v 2 

has "been chosen so that the leading term ^ — ty + 

I 2 \ 2t 

satisfies the condition of equation (124b). The constant 
k x is then determined so that. the coefficient of. S^Cx) 
in equation (124b) vanishes and 




61 



{(I s + 3I W ) 4- ^£ (i + 3i )'v (127) 
v. 2tw J 



after which the constants k 2 , 
mined so that the coefficients 



k 3 , can he deter — 



of M I7 (x), M VJ -( X ), 



in equation (124b) vanish. It appears that, if H(z) in- 
volves only odd powers of x, the polynomial P(x,y) 
satisfies also equation (C2) as well as equations (Cl) , 
(C4) , and (C5). It may also he noticed that the leading 

M(x) h / A-A 
term . ... .... — ( y + — — ) corresponds to the stress distri— 

I 2 \^ 2 1 / 

hution given by the elementary theory of the strength of 
materials, according to which the sp*anwis e nor mal stress 
a r does not vary in the transverse direction. 

The conditions of equations (123c) and (123d) are now 
fulfilled if Y n (y) satisfies the Sturm-Li ouville^ problem 
consisting of the differential equation 



w 



Y„ = 0 



n n n 

and the homogeneous boundary conditions 



(128) 



Y (w) + w -S Y n '(w) = 



r 



(129) 



If it is assumed that f (y) has a continuous second 
derivative in (o,w), it follows from the remarks made 
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in connection with equations (100) and (101) that, if 
condition (97) is imposed., equations (123a) and (123b) 
imply 

x n<°> = - f *(o,y)Y n (y) dy (l30a) 
o 

lQ X n (0 + ax X a ' (D = (y)-a 0 P(l 1 y)~a 1 P 3c (l,y)|y tl (y)dy (130b) 

o 

and these conditions, together with equations (122a) and 
(122c), determine X n (x). 

Hence, finally, the SQ-lution of the problem of equa- 
tion (c), in case M(x) can he oxpressed as a polynomial, 
is given by equation (120), whore P(x,y) is given in 
equation (125), Y „(y) is determined "by equations (128) 
and (129), and X n (x) is determined by equations (122a) 
and (130). 

The "effect ivo area" of the eorabined cover sheet — and 
central stiff ener, defined by the r elati onship 



A eff a xf x »*3 = ( A m + 2tw) eff .a x (x,w) = 2t J* cr^dy + A Q C7 2 (x , o ) (131) 

o 



can be expressed in terms of the elementary normal stress 

and the actual flange and stiffener stresses as follows. 
If equations (7a), (88), (90), and (126) are used 

A Rff ct 2 .(x,w) = 2t [H(x,w) - H(x, o+)] + A n ct x (x,o) 

= 2tw ^ ~ a b (x) - 0 a x (x,w) J 

= 2tw [O + 1) a h (x) - B a x (x,w)] 

and, finally, 
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L eff 



r 

= 2tw < 1 - O + 1) 



I 



1 - 



r 



(132) 



In particular, if no central stiff enfer is present, the 
ratio of the effective width of the cover sheet to the 
actual width is given by the expression 



w 



1 - 



a (x,w) 



(133) 



The evaluation of the solution may he illustrated 
by taking "the case of a concentrated tip load 3f , f or 
which - 



M(x) = Ex 



(134) 



and the stress function is ■ determined for an end— support 
condition of the form 



(135) 



It is further assumed ±hat no. central sheet stiff oner is 
prejont. ffirst, from equation (125), it follows that ' 



?U,y) = cr Q f y 



CT o = 



Phi 
31 



(136) 



Next, from equations (128) and (129), the constants \ s 
are determined as the solutions of the equation 

tan X n + p\ n = 0 " 



(137) 



while the functions Y n (y) take -the form 
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Y n (y) = V n Bin \ n L 

Y 3 = - _ 2 - - 
° w(l+ B OOfi 8 X n ) 



(138) 



The terminal conditions governing the functions X n (x) 
tocoae , according to equation (130), 



X n (o) 



± 0 X n {l) + a. i X n '(l) = ~0 a 0 +a 1 )a o Y n 0+ l) 



v > (139) 
cos X„ f 



J 



p„nd hfincfi, from equation (l 22 ) 

cos X n < l a 0 -+ai)sinh u-^ ^ 

X n (x) = v 2 C7 0 Y n (6+l) 



^n I a 0 sinh p, n +ai |J. n cosh p. n 



(140) 



If these expressions are introduced into equation (120), 
the stress function is 'deter ain-ed in the form 



H = w ct , 



5 Z + 2(8+l>> 



21, cosX n sin\ n ^ (I a 0 + a 1 )sinh u- n y 



r 



-> \ n (l+Bcos 3 \ n ) la 0 sinh(i n + ai |A n coshn n 
n= l 



(141) 



The stress patterns derived from equation (141) in the 
cases a 0 = 0 and- a x = 0 were discussed in a preceding 
sect ion. ' 

In case the sheet, the stiff ener, and the side struc- 
ture are rigidly clamped at x = I , so that a 0 = f (y) 3 0 
in equations (C3) and (123a), the general solution for arbi- 
trary M(x) can also he obtained . very 3imply by the alter- 
nate procedure mentioned at the beginning of this section. 
If P(x,y) s 0, equations (l23a)'and (123b) are satisfied 
if X n (°) = X u '(l) = 0 so that, from equation (122) , 



^n = 



(2n - l)iTi 



X n (x) = sin 



(2n - ;1)ttx 
21 



> 



(142) 
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Equations (122) and (123c) require that 



0 



T B (7) = c n ^sinh 8 n t + ---- - n 

(2n - 1)tt /i w 
2 V 0 J 



8„ cosh 5 



8 n = 



(143) 



and , if 



M(x) - = / a tt sin ~ 



(144) 



n= l 



equation (123d) determines H in the form 
= S -- — ( sinh6 n - + — — cosh 5 n - j 



H 

where 



sin !^Zi2^_(l45) 
2 I 



A 



^n = I s (^sinh 6n + 8 n cosh 6 n 



+ 8 n I w (c 



S n ^ (145a! 



cosh 6 n + — r~ 6 n sinh 5 
2tw 

This solution was given for the special case in which. A. m = 0 

and M = M 0 sin ~ "by Younger (reference 4) and for arhi- 
v 2 1 

trary M with A m = 0 hy Kuhn (reference 3). 

The Semi—Infinite Sheet with One Stiff ener 

In order to determine the stress function correspond- 
ing to the state of stress in a s emi— inf ini t e sheet sub- 
jected to a load introduced into a stiff ener normal to the 
straight-edge (fig. 3(a)), a different type of procedure 
is needed. If the y axis is taken along the straight- 
edge and the X axis coincides with the stiffener and if 
the sheet is assumed rigid in the y direction, the stress 
function H is determined "by the differential equation 
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ox s B oy 3 

and the . boundary conditi-ons 

H(o,y) = 0 (147a) 

H(»,y) = 0 (147D) 

_2t H(x.o) + A = P (147c) 

by 

H(x,«>) = 0 (1474) 

If it is assumed that H is the Laplace transform of a 
function g to "be determined, it follows that 

H(x,y) = J e~ yP g(x.P') *P < 148 > 
■ o 

The differential equation (146) requires that g satisfy 
the equation : - 



ox 2 



£s + k 2 p 3 g = ;o 

•J E 



(149) 



while condition (147a) is satisfied if 

g(o,p) = Q (150) 

so that 

g(x,p) = h(p) sin; k x P (151 ) 

where h is an arbitrary function of p. 

If this expression is introduced into equation (148), 
'condition. (147c) requires that h(p) satisfy the integral 
equation - 

oo 

/ (2t + A p) h(p) sin k x p dp = -3P (152) 
which has as solution the expression 
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h (p) = _ 2T ■ _J. _ (153) 

17 P (2t + Ap) 

Thus a formal solution of the problem is determined 
in the form 



CO 



E =, - H [ 6 ~™ J^2JL£P dp ( 154 ) 
ttA / / ^ 2tA 



It can be shown without difficulty that the conditions of 
equations (l47b) and (l4?d)'are satisfied and, furthermore! 
that the expression (154) constitutes the rigorous solu- 
tion of the problem. 

The stress distribution associated with equation 
(154) was discussed in a previous section, .. . 

COMPARISON OP APPROXIMATE PROCEDURES 
Axially Loaded Panels 



Although the analysis of panels is not included in 
reference 1, approximate solutions for most of the cases 
considered in the third section can be found by the pro- 
cedure of reference 2. In that procedure the effect of 
the transverse stretching of the sheet is neglected, in 
accordance with the present theory. In addition, it is 
assumed that the resistivity of the sheet in the spanwise 
direction can be neglect-ed if the sheet area ia added to 
the area of the central longitudinal, or treated as a 
fictitious central stiffness in the absence of a longi- 
tudinal, and a substitute width of one— half the actual 
panel width is used in the subsequent calculations. 

In order to investigate the agreement between the 
results of this procedure and the present results, panel 
8 (fig. 4(h)) is considered. In this case the expres- 
sions' for the flange stress a x (x,w) and for the average 
spanwise normal stress cr a (x). are obtained in the form 

2P / A m + 2t w _e x\ . . 

cr x (x,w) = — (l + — e K w ) (155) 

A m + 2A 6 + 2t w V 2A e / 

an d 
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«r a (x) = JSL — (l - e K W ) (156) 

A m + 2A e + 2t w V / 

where K is a alimens ionles s constant defined by the 
relationship 

K 3 = — t —f 1 + JS\ (157) 

B 2 V + 2t v A e / 



The transverse variation of the spanwise normal stress 
at each sect ion. is assumed in reference' 2 to follow a 
hyperbolic— cos ine law, in terras of whiob the t'.tjrees vari- 
ation in the longitudinal can "be expressed. 

"8 or the numerical evaluation of the solutions two 
cases are chosen ... 



(a) 



00 



Am 
A_ 



= t w 

1 . 

= — t w 

2 



.-I 



A m = 5 t w 



A £ 



(158) 



(159) 



a = 



in- 



corre sponding , 
5 
2 

present a 
It may he 
estimates 
est imates 
crease in 
smallef as 



respectively , 
the third section. 



to the cases 
Figures 16(a) 



a = 



1 
5 



and 
and 16 (b ) 



comparison of the results of the two procedures, 
noticed that the Euhn solution in general over— 
the magnitude of the flange stresses and under— 
the longitudinal stresses, with a resultant in— 
the predicted shear— lag effect , which becomes 
the contribution of the sheet to the tot~al area 



of the panel . decreases . * 

The deviation of the Kuhn solution from the exact 
solution in the neighborhood of the loaded end of the 
panel is partly due to the fact that, in this region, the 
transverse variation of the normal stress cannot be sat- 
isfactorily approximated by a hyperb olic— c os ine law. 
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Doubly Symmetrical Box Beam with Concentrated Tip Load 



If' it is assumed that no central .sheet stiffener is 
present, the stresses given "by the procedure in reference 
1 can be .written in the form 

~> 



a-Jx.o) = a, 



1 (3E + D(6B+ 1) sinh Kr T 

2 51p 2 + 18 8 + 2 K E cosh E H 



(160) 



r 



£ + 7 



6.3 + 1 



s inh K- 



•R 



51 6 s + 18 8 + 2 Kr cosh Kr 



(161) 



where a 0 and 8 are defined by equations (44) and (45) 
and the parameter K E is defined by 



21 



G (P + 1) (6 8 + 1) 
3 51 8 S + 18 8 + 2 



0) 



(162) 



This solution was obtained by an application of the method 
of least work, under the assumpt i on of a parabolic trans- 
verse distribution of the spanwise normal stress ^x* 



According to the procedure of reference 2, it is nec- 
essary t o .incorporate the sheet area .2t w into ^ ficti- 
tious central stiffness. The foilowing expressions are 
then obtained for the flange stress and the average sheet 
stresst • . _ 



o- x (x,w) 



r 



1 sinh K 7 



\ 



8 K cosh -K j 



(163) 



a a (x) 



= a . 



1 



s inh K 



K cosh K 



fl 

r 



(164) 



wher e 



(165) 



The stress variation along the center line of the sheet is 
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again determined "by assuming that the sheet normal stress 
varies according to a hyperb olic— cosine law in the trans- 
,versa direction. 

lor a numerical comparison of" these solutions, the 
following values of the parameters are chosen in the same 
way as in the third section: 



r: =» 5 



0 = -r- = 1 



(166) 



G- 



3 
8 



The shear— lag 

compared with the 
present procedure 
Heissner solut-ion 
solution, except 



functions 



( 

according 



J 



(x,w) 

o-o 
to the 



-f) 



and 



two procedures , are 



the 



corresponding functions obtained "by 
in figure 17. .It is seen that the 
is in close agreement wit-h the present 
in the immediate vicinity of the fixed 
end, where differences of 11 percent and 8 percent occur 
in the stress corrections at the center of the shoot and 
at the flanges, respectively, corresponding to differences 
of 3,7 percent and 0,8 percent, respectively, in the val- 
ues of the stresses themselves." The shear-lag eff~e~ct 
predicted by the Kuhn solution is again considerably 
greater than that given by the present- procedure, 
ences of 62 percent and 35 percent are present in 
root-stress 'corrections at the center of the sheet and at 
the flanges, respectively, corresponding to differences 
15 percent and 4 percent in the actual stress 



Differ- 
the 



of 



values 



In view of the fact that for a uniform box beam sub- 
jected to a concentrated tip load, the shear— lag effect 
is comparatively small and is appreciable over only a 
small part of the span., it seems desirable further to com- 
pare the .solutions obtained by the approximate methods and 
by the present procedure for other types of beam and load- 
ing where the shear— lag effects are more pronounced. 



Effect of a Out— Out in the .Cover Sheet of a Bos Beam 



The Heissner procedure is not directly applicable to 
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the analysis of c over sheets with cut— outs , since in. the 
vicinity of the cut— out the transverse distribution of 
cr x is no longer approximated "by a paraoola. It may "be 
pointed out that a least— work solution of the type con- 
sidered in that procedure would, "be obtained if it were 
assumed that the difference "between, the actual normal 
stress o* x and the elementary stress o* 0 ^ is expres- 
sible as the product of a function of x and a function 
of y. If then the conditions of equili"br ium were satis- 
fied and the method of least work were applied, a correc- 
tion would he 'ohtained that would be identical with the 
first term of the series given in equation (57). The cor- 
responding stress pattern, for the values of the parameters 
listed in equation (166), is compared in figure 18 with 
the pattern obtained "by retaining the complete series. 
It is seen that the first term of the series affords a 
reasonably accurate correction at a distance from the ' 
cut— out greater than one— half the width of the sheet. 

The corresponding solution according to the Euhn .pro- 
cedure can "be written in the form 



f sinh E f 



sinh E 




s'inh E f- 



sinh K 



(167) 



(168) 



where,, as before, o" a (x) is the average sheet stress and 
-E is defined "by equation (165). The stress patterns", 
for the values of the- parameters given in equation (166), 

are compared, with the present results, in figure 19* It 

is seen that the shear— lag effect is greatly overestimated 
by the Euhn solution and that, except in the immediate 
•vicinity of the cut— out where the stresses are already 
known, a hotter approximation to the actual stresses is 
given by the first two terms of the present solution. 



NUMERICAL APPLICATION OT THE THEORY TO A TEST PANEL 

In order to illustrate further the numerical evalu- 
ation of the exact theory, the panel sketched in figure 
4(j) is analyzed. The results are then compared with the 
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experimental test results reported in reference 3 and with 
the approximate results pr e sont 3 4 in that" paper. 

If the stress function (30) is modified in accordance 
with the st-a-tement at the beginning of- the t-hird soction 
and if : 



and 



»a = 0 



a = 2a 



ra 



the expressions, for t ho . s t if f e nor and flango strossos fol- 
low from equation (7a) in tho form 



ct (x,o) = 



CT x (x,w) = 



cosh/ ^ X n -b 3 
v 3 u w 



(3a m + l)tw 



'■ — 1 K 



cosh '/— X„ — 



(169) 



• o 



(3am + l)tw 



as 



1 *r 



K 



f_ Tl 7 n ( cos K n ~ a m x n 8 x n) — 
n= 1 cosh 



cosh /— X 



l-x 



>(170) 



whore, from equation (34), the constants X n arc tho 
solutions of trho equation '! ' J 



tan X n = 



3ajE X c 



(171) 



and tho definitions of tho remaining constants, as given 
in equations (32) and (33), tako the form 



(2X n - sin 2X n ) - 2a m X n (cos 2X n - 1 ) 

+ a m S X n 3 ( 2X n + s 2X n )V (172) 

-\ 
j 



Vn 3 = 4X n 3 \ 
^n~ w r 

F£ = 7T' l( sin X n - X n cos X n ) + a ra (l + a ra )X n sin X^ (l73> 
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In accordance with, the- notation of t he.. .third sec- 
tion, the following data are to he assumed: 



■* m 




0 








w 




4.60 


I 




36 


t 




0.016 


A e 




0.366 


G 
S 




0.4 



and hence, from equation (31), 

a = 2a. 



n 



a a = 2.487 



(174) 



(175) 



First, fron equation (171) , the parameters X are deter— 



rained readily by a method of suc'cessive appr oxinat ions , 
after which the following table is ' cons truct ed fron equa- 
tions (172) and (173) : 



n 


*n • 


K n v 3 


cos K n - a n \ n sin \ n 


1 


0.7521 


1. 0522 


-0.547 8 


2 


3.323 


-.1229 . 


.5704 


t-r 

o 


6.397 


.0412 


-. 8134 



(176) 



The results obtained by introducing the data of aquations 
(174), (175), and (176) into equations (169) and (170) are 
presented in tabular f orn and the computed stresses are 
compared with the corresponding values determined by 
approximate formulas givon by Kuhn in reference 3 , as fol- 
lows : 
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x 

0 
3 
5 
8 
10 
15 
20 
25 
30 
36 

x 

0 
3 
5 
8 
10 
15 
20 
25 
30 
36 
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Solution "by equation (169) 



First 
tei*m 



3850 
3850 
3850 
3850 
3850 
3850 
3850 
3850' 
3 850 
3850 



Sec ond 
term 



-4050 
-2970 
-2420 
-1770 
-1440 
-860 
-530 
-330 
-2,30 
-190 



Third 
t or m. 



4-70 
120 
50 
10 



Fourth 
term 



-160 
-10 



cr x (x,o) 



110 (0) 

990 
1480 
2090 
2420 
2990 
3320 
3520 
3620 
3660 



Kuhn 
s olut i on 



a x (x,o) 



0 

990 
1500 
2110 
2420 
2960 
3300 
3500 
3600 
363.0 



Solution by equation (170) 



Kuhn 
s olut ion 



First 
term 



3850 
3850 
3850 
3850 
3850 
3850 
3850 
3850 
3850 
3850 



S econd 
term 



2220 
163 0 
1330 
970 
800 
480 
290 
190 
130 
110 



Third 
term 



270 
70 
30 
10 



F ourth 
term 



.130 
10 



tf x (x,w) 



6 4 7 0(6560) 

5560 

5210 

4830 

4650 

4330 

4140 

4040 

3 980 

3960 



cr x (x ,w) 



5960 
5420 
5140 
4810 
4640 
4330 
4160 
4050 
4000 
3980 



(177) 



(178) 
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The calculated stresses, t ogether with the test data 
taken from reference 3, are presented graphically in fig- 
ure 30. 

It is seen that the convergence of the series if. 
sufficiently rapid that only two terms of the correction 
series need he retained. While the prescribed "boundary 
values (noted in parentheses in equations (177) and (178)) 
are not yet attained, this fact is of no consequence since 
the values are known and since the calculated stress val- 
ues at a distance of only 3 inches from the free end are 
apparently accurate to three significant figures. The 
Kuhn solution, was calculated "by associating one— half the" 
sheet area with the longitudinal and the remainder with 
the two flanges, and is in very good agreement wTth the 
present solution except near the point of load applica- 
tion. 

The discrepancies hetween the theoretical and the 
test values, as was pointed out in reference 3, can be 
attributed to two factors: (l) because of the magnitude 
of the loading, an observed tension fold was developed on 
each side of the panel near the loaded end, reducing the 
effective shear modulus in this region; and (3) the fact 
that the flange stress readings near the root of the panel 
are actually loss than the statically requisite limiting 
values for an infinitely long panel must be accounted for 
by assuming an elastic deformation of the steel triangle 
which served as the panel support. 



MISCELLANEOUS EXTENSIONS AND MOD IE IOAT IONS OP THE THEORY 



The procedure given in this paper can be extended in 
many cases to apply to box beams or panels the cross- 
sectional characteristics of which vary along the span. 
Three explicit cases may be mentioned here. 

Tap^r in baam h nig ht.— The basic equations (0) again 
apply. If it is assumed that the side— web characteristics 
also vary in such a way that the ratio of the contributions 
of the web and of the sheet to the total stiffness of the 
beam remains constant along the span, so that 



Analysis of Tapered Box Beams 



I 



w 



(179) 



3 = 



I 



= constant 



s 



56 NAG A Technical Note Bo. 894 



tho effect of varying height can he taken into account by 

cL n / H (x) h" 



u j 21 

merely replacing. — M(x) 

21 dx 11 



■by 



in tho 



definition of the function P(x,y) as given in equation 
(125). 

Tavar in cover— sh e et thickness.— If the problem is 
formulated in terms of the stress rbsult"ants rather than 
in terms of the stresses themselves and if a stress func- 
tion 3> is defined, in terms of which 



txr, = _ 

t T .-|» 
OX 



(1,80) 



it is found that 
equat i on 



$ is determined. by the differential 

d 2 $ 



bx It ox J B 



(181) 



and "boundary conditions that differ from those which 
apply-to the function H in equation (C) only by the sub- 
stitution of i M(x) for M(x). If it is assumed that 



A,. 
I 



const ant 



w 



— = constant 

I a 



(182) 



the solution can be readily obtainod in terms of— Bessel 
functions t by a modification of the procedure of a preced- 
ing sect-ion, for beams in which the cover— sheet thickness 
varies according to a lav/ of the form 



t (x) = t Q (x + x Q ) 



(183) 



For example, if a box beam- in which 
thickness trgpors linearly to zero 'at the 
that * 



the cover 
free end , 



—sheet 
s o 
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(184) 



is subjected to a uniformly distributed load of intensity 
v _ , so that 



M(x) = | p 0 x 2 



(185) 



and if the "beam is fixed at the end x ~- I (no central 
sheet stiff ener "being present), the stross function $ 
is obtained in the form 



+ .4 0 + 



i) /a ^ 

* Gr I 



cos \„ s 
n 



in X a £ Ix^lfcn f 



n= i 



\ n s (l+0 cos 3 X a ) I o(l J 'n) 




(186) 



and Ij are modified Beasol functions of the 



where I 

first kind, of orders zero and one, o* 0 is the maximum 
elementary , stress , and the remaining quantities are de- 
fined in equations (45) and (46). S'roa equation (180) 
the stresses are given "by the expressicTns 



CT x = CT o 



I 



4(6 + 1) 




7 

Ma v 



cos \ n cos \„ — IV — 



Gr I 



> 

n^i M 1 * 0 cos s X n ) 



W 



£+4 (0+1) y 

w i 



w 




(187) 



(188) 



An approximate solution of the same problom is given in 
reference 5. 



Tapered cover— sheet width r — She principal modifica- 
tion necessary in case the width of the shoot varies in 
the spanwise direction arises from the fact that the nor- 
mal stress in the flange, xjan np longer be identified with 
the spanwise normal stress in the adjacent sheet fibers, 
"but is related to the sheet stresses "by the expression 
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[ff(x)]*-, = o- y (x,w)cos 3 6 + cr (x,v)9in s 5 +T(x,w)sin 26 (189) 

i lange . ■ - * 

where 8 is the angle "between the x— axis and tho normal 
to the edge of the sheet. 

A More Ooneral Class of End-?Suppor t Conditions 

If equation (03) is replaced "by a condition of the 
more general form 

a 0 H(i, y ) +ai Mil^i +a3 + a 3 5!5iV^ + ....£( J ) (190) 

ox o.-t 3 ox 

equat-ian (122) can "bo usod to rewrito tho equation replac- 
ing equation (l23"o) in the form 

V{^a 0+ a 3 (^) 2 + ]x»(0 

+ [a, + a 3 (^) S + .....]x n «(l)| Y Q (y) 

= f(y)- Ca 0 P(l,y)+a 1 P x (J,y) + a s P. x . x (i,y) + a 3 P xxx (^3r)+...](l9l) 

The only modification in the solution therefore consists 

in'replacing the terminal condition of equation (130b) by 
■the condition 

+ a3 (^)%....]x n (0 + [a 1 + + ] X n'( J ) 

= / |f 00 - a 0 .PO','u)-a 1 P I (i l u).-a. P xx< l ' u >- }? n (u)du 
° (192) 

Jor example, if the "box hear* is attached to an elas- 
tic springlike support, the end condition is of the form 

• u (1.7) =— c.cr x (l,y) ' • ■ (193) 

If both sides of equation (193) are" differentiated with 
respect to y, and if equations (5), (.7), and (Cl) are 
used, this condition can be written 
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ox 



ox 2 



= 0 



(194) 



lor the case of a concentrated tip load, equation (139b) 
is now r eplaced , according to equation (192), by 

cB (Ty X n(0 + X n «(0 = — cr 0 Y n (P + 1) ~ x — * (195) 



so that, in place of eauation (140), 

cos X„ 



X n (x)=w 3 a 0 Y n (6+l) 



smn M-n y 



*n c E s 

u- n sinh (a, n + u- n cosh n n 



(196) 



and the stress function is ob.tainod in the form 



H . . ,„ L 



+ 3(0 + 1) A f > 



cob X-rt ain- X-n — 



w 



.sinh n n j 



G 1 rfr'i ^ (1+P cob3 \i> fi-l^sinhHtf-cokhp.: 



> (197) 



Box Beam with One Cover Sheet . 

If a box beam carries a cover sheet on only one side, 
the equation of moment equilibrium about the neutral axis 
becomes (see reference 6, p. 11, equation (12)) 



I w a z (x,w) + 2t(d3 + h. ) J' a s dy+ A m (d 



3 + M(x)d (198) 

A w y 



where, in addition to the cons'tants previously defined7 
A w is the cross-sectional area of the two side webs 
(including the corner flanges) and d is the distance from 
the plane of the cover sheet to the neutral axis. If epila- 
tion (198) is expressed in terms of H and if equation '_ " 
(88) is taken into account, epilation (90) is now .replaced 
by the condition \ 
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It follows that the results of the preceding sections can 
"be applied directly to this case if I s is replaced by 



The shear— lag analysis of the present paper presents 
mathematically exact solutions of several problems that 
have been previously treated in general by approximate 
mothods. Although the rosulte .strictly apply only to 
beams and. panels that are theoretically rigid in the 
transverse direction, it seems probable that they are 
applicable with reasonable accuracy in actual cases when 
stiff chordwisc ribs are presont. The solutions aro ob- 
tained in the form of rapidly convergent infinite series 
that are much noro adaptable to numerical computation 
than exact solutions that have been given elsewhere in 
certain special cases » 

On the basis of comparisons of— the present results 
with approximate solutions given by Reissner and by Kuhn 
and Chiarito, it 'appears that--in cases when both the 
approximate solutions are applicable, for example, uniform 
rectangular box beams without cut— outs, the Reissner solu- 
tions are in better agroement with the exact solutions than 
aro the solutions given by the methods of Kuhn and Chiarito, 
In other cases for which the Reissner procedure was not 
designed, for example, box beams with cut— outs and panels 
leaded by concentrated axial -f or ces , the solutions given 
by the procedure of Kuhn and Chiarito predict shear— lag 
effects that are in general considerably larger than those 
given by the present procedure. 
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Table 1.- Solutions of tan ^ + oO> n « 0, 



n 




> 


n 




<X= 1 


dm 5 


1 ' 


2.0288 






1.6887 


2 


4.9132 






4.7644 


3 


7.9787 






7.8794 


4 


11.0866 






11.0157 


8 


14.2074 






14.1613 


6 


17.3364 






17.2903 


7 


20.4692 






20.4301 


a 


23.6043 






23.6704 


9 


26.7409 






26.7110 


• 10 


29.8786 






29.8618 


u 


33.0170 






32,9928 


12 


36.1660 






36.1339 


13 


39.2064 






39.2750 


14 


42.4361 






42.4162 


15 


46.6760 






46.6676 



Table' 2. - 8treef*a in Panel So. 1. 



X 

* 


CI- 1 




a = 6 


a - 


0 




















0 
0.1 - 
0.2 
0.4 
0.8 
1.2 
1.6 
2.0 
5.0 
5.0 


0 
.036 
.071 
.139 
.266 
.343 
.402 
.440 
.482 
.499 


i. 

.86 

.816 

.728 

.688 

.676 

.546 

.628 

.508 

.601 


-.79 

-.643 

-.338 

-.172 

-.097 

-.068 

-.034 

—010 

-.001 


0 
.062 
.113 
.213 
.391 
.330 
.628 
.697 
.784 
.828 


1. 
.97 
•951 
.926 
.891 
.871 
.868 
.860 
.839 
•834 


-1.06 
-.806 
-.653 
-.319 
-.200 
-.129 
-.085 
-.030 
-.004 


0 

.096 
,190 
.367 
.647 
.819 
.912 
.968 
.994 
1.000 


10.4S' 

5.26 
2.726 
1.546 
. 1.221 
1.097 
1.044 
1.006 
1.000 



I 



Table 3.- Solutions of tan An =• — * P - r 

c&a 3 + 1 



n 




ff= 0 


a= 1 


a~ 5 


1 


4.4934 


3.4056 


3.2028 


2 


7.7253 


6.4538 


6.3147 


3 


10.9041 


-—9.5272 


9.4469 


4 


14.0662 


12.6448 


12.6823 • 


6 


17,2208 


' — 15,7710 


15.7007 


6 


30.3713 


18.9023 


18.8602 


7 


23.5195 


22.0364 


82.0002 


8 


26.6661 


25,1724 


26.1407 


9 


29.8116 


28.3096 


28.2814 


10 


32.9564 


31.4477 


31.4223 


U 


36.1006 


34.5864 


34.6633 


12 


39.2444 


37.7266 


37.7044 


13 


42.3879 


40.6632 


40.8456 


14 


46.5311 


44.0060 


43.9369 


15 


48.6742 


47.1461 


. 47,1281 



? 

o 

.... . . » 

Table 4. - Stresses In Panel Ho. 2. g 



X 
V 


a« 1 


a- 6 


a= 0 


fax?) 


|rU,o) 




tew 




fTU.w) 




^T(X,0) 


0 

0.1 
0.2 
0.4 
0.8 
1.2 
1.6 
2.0 
5.0 
5.0 


1. 

,911 
.867 
.818 
.776 
.761 
.768 
.753 
.751 
.750 


.21 

,189 

.160 

.113 

.055 

.026 

.011 

.005 

.001 

.000 


-.55 

-.340 

-.171 

-.039 

-.024 

-.010 

-.004 

-.001 

-.000 


1. 

.980 
.972 
.959 
.944 
.941 
.939 
.937 
.936 
.936 


.259 
.224 
.192 
.140 
,070 
,054 
.016 
.007 

.001 

.000 


— Oft 

-.67 

-.420 

-.326 

-.065 

-.036 

-.016 

-.007 

-.001 

-.000 


11.6 
6. 57 
4.225 
5.279 
5.081 
5.026 
3.. 006 
3.001 
3.000 


m 

l.oie 

.722 
.416 
.150 
.055 
.018 
.006 
.000 
.000 



Table 5.- Solutlona of oot Xn - So&a * 0. 







ft 




tt« 1/2 


a- 5/2 


1 

X 


.86053 


.43884 


2 


3.4206 


3.8039 


3 


6.4573 


6.3148 


i 


9.6295 


9,4469 


S 


1J2.84SS 


12.6883 


6 


16.7718 


16.7807 


7 


18.9084 


18.8608 


8 


22.0366 


88.0002 


9 


25.17B4 


86.1407 


10 


£8.3096 


88.8814 


11 


31.4477 


31.4883 


18 


34.6664 


34.3633 


13 


37.7866 


37.7044 


14 


40.8668 


40.8466 


15 


44.0060 


43.9869 



labia 6.- Stress* a In Panel Ho. 6* 



1/3 



a= 5/2 



0 

.036 
.071 
.106 
.139 
.866 
.343 
.408 
.440 
.488 
.499 



1. 

.88 

.616 

.766 

.788 

.688 

.676 

.646 

.688 

.808 

.601 



-.79 

-.643 

-.419 

-.338 

-.178 

-.097 

..068 

-.054 

-.010 

-.001 



0 
.068 
.113 
.164 
.813 
.391 
.630 
.688 
.697 
.784 
.888 



1. 

.97 

.961 

.937 

.983 

.891 

.871 

.868 

.860 

.859 

.834 



labia 7. - StreaBOB In Penal Ho. 7. 



x_ 


a. 


•= 1/2 


CC: 


• 5/3 


>*Ox(x,») 






I s ltx,») 


0 


1. 


*■ mm 


1. 




0.06 


.867 


..964 


.966 


-1.86 


0.1 


.786 


-.718 


.980 


-1.606 


0.8 


.667 


-.600 


.868 


-1.347 


0.4 


.608 


-.319 


.766 


-1.091 


0.6 


.400 


-.831 


.683 


-.940 


0.8 


.319 


-.177 


.614 


-.888 


1.0 


.867 


-.159 


.668 


-.738 


1.6 


.161 


-.080 


.483 


-.361 


8.6 


•063 


-•028 


.849 


-.330 


6.0 


.004 


-.008 


•006 


-.006 



■ 

o 

p 

H 

s 



Tabl* 8a. -Stresses In Panel Ho* 8. (<-£) 0 



X 

w 




K<4> 








0 


0 


0 


1. 


.09 




0.06 


.008 


.018 


.876 


.09 


-.869 


0.1 


.016 


.053 


.800 


.087 


-.689 


0.16 


.082 


.063 


.744 


.084 


-♦603 


0.8 


.030 


•069 


.697 


.082 


-.419 


0.3 


•046 


•101 


.625 


.077 


-.312 


0.4 


.061 


.188 


.669 


• 074 


-.845 


0.8 


.088 
•113 


.178 


.488 


.067 


-.164 


0.8 


.801 


.432 


.060 


-.117 


1.0 


.336 


.220 


.392 


.068 


-.087 
-.083 


2.0 


.807 


.848 


•806 


.088 


5.0 


.848 


.260 


.862 


.001 


-.001 



i 



Cable 6b. ■ 



Streeeea In Penal Xo* 9. (ot» 5] 



0 

o.oe 
o»i 

0O5 

0.2 

0.3 

0,4 

0.6 

o.e 
l.o 

CO 
5.0 



0 
.007 
.010 
.024 
,082 
.048 
.062 
.093 
.122 
.149 
.256 
.383 



0 
.031 
.060 
.082 
.100 
.101 
.196 
.206 
.314 
.349 
.408 
.416 



1. 

.963 

.936 

.913 

.893 

.858. 

.829 

.778 

►736 

.701 

.580 

.450 



.40 
.40 
.400 

.397 
.392 
.386 
.370 
.350 
.326 
.210 
.043 



-1.46 
-1.200 
-1.066 
-.960 
-.805 
-.705 
-.570 
-.479 
-.411 
-.221 
-.044 



Table 9, -Wga S treason in UnatlfTeneiL Panela. (tt. q) 



w 

~3 — 
oa 

0,2 

0.3 

0.4 

0.6 

0.8 

1.0 

1.25 

1.00 

2.0 

2.0 



Panel No. 7 



10.524 
6.452 
3.842 
5.092 
2.442 
2.194 
2.068 
8.032 
2.002 
2.002 
2.000 



Panel So. 8 



10.330 
5.072 
3.280 
2.308 
1.400 
.900 
.596 
.364 
.224 
.086 
.032 



Panel No. 2 



11.612 
6.572 
4.973 
4.223 
3.564 
S.279 
3.148 
3-071 
3.034 
3.008 
3.002 



Eablo lO.-Streesea In Carer Sheet of Box Beam Without Cutout. 



X 

T 


P 


= 1 


P ■ o 










0 


0 


0 


0 


0 


0,1 


.099 


.ioo 


.100 


.100 


0.2 


.199 


.201 


.200 


.200 


0.3 


.297 


.301 


.300 


.300 


0.4 


.594 


.403 


.400 


.401 


0.5 


.439 


.506 


.499 


.502 


0*6 


.580 


.609 


,696 


.606 


0.7 


.054 


.710 


.089 


.712 


0.8 


.734 


.831 


.772 


.833 


0.9 


.782 


.959 


.833 


1.000 


0.95 


.795 


1.033 


.860 


1.150 


0.97 








1-.258 


0.99 








1.487 


1.00 


.800 


1.120 


.868 


CU3 



o 

p 

H 

Table U.- Eranaverse Diatrihutlon of Hoot Bonwl Streas as 
In Cover Sheet 'of Box Bean Without Cutout. ° 





. Oxd.yVoo 


P - 1 


P - 0 


0 


.800 


.850 


0.1 


.801 


.856 


0.2 


.807 


.866 


0.3 


.817 


,880 


0.4 


.831 


.900 


0.5 


.851 


.928 


0.6 


.876 


.966 


0.7 


.908 


1.080 


0.8 


.951 


1.100 


0.9 


1.011 


1.242 


0.99 




1.720 


0.999 




2.198 


1.000 


1.12o" 





Table 12, - Spanwlse Variation of Cover-aneet Effective Mdtk. 



Fixed sheet 





P - 1 


p « 0 


P - 1. 


p " 0 


0 


1. 


1. 

1.000 


1. 


1. 


o.x 


.994 
.993 


.987 


.997 


o.a 


1.000 


.971 




0.3 


.991 


•990 


.930 


«9QS 


0.4 


.988 


.998 


.926 


.986 


0.6 


.981 


#0Qd 


.889 


• 96B 


0.6 


.971 


.998 


.633 


.932 


0.7 


.934 


.983 


.743 


.866 


O.B 


.9126 


.960 


.606 


.701 


0.9 


.877 


.900 
.826 


.393 


.421 


' 0.96 


.839 


.248 


.227 


0.97 




.771 




M 


0.99 




.666 


IP 




1.00 


.786 


1 0 


0 


0 



Free sheet 



labia lS.-Streasea In Cover Sheet of Box Beam 91th Cutout. 



X 

\ 


P ■ 


i . ! 


P - o 










0 


0 


0 


0 


0 


0.1 


.096 


.101 


.100 


.100 


0.2 


.191 


.204 


.199 


.201 


0.3 


.281 


.308 


.298 


.302 


0.4 


.365 


.416 


.394 


.406 


0.6 


.454 


,629 


.484 


.616 


0.6 


.479 


.666 


.668 


.644 


0.7 


.479 


.803 


.694 


.818 


0.8 


.408 


.997 


.646 


1.142 


0,9 


.842 


1.292 


.347 


2.137 


0.96 


.127 


1.53 


.186 


4.188 


1.00 


0 


2.000 


0 


M 



i 



I 



Table 14. 



Stresses in stiffened 


Stresses in unntiff enad 




half-sheet , 


half-sheet'. 








Vft 




0 


1. 




0 




0.01 


.968 


8,674 


0.01 


63.66 


0.10 


.882 


1.188 


- 0.10 


6.566 


0.26 


.685 


.717 


0.85 


2.547 


0.60 


.348 


.428 


0*60 


1.273 


1. 


.396 


.219 


1. 


.637 


2. 


.264 


.098 


2. 


.318 


3. 


.120 


.022 


6. 


.127 


10. 


.063 


.006 


10. 


.064 



i 

e 

o 



labia 15.- CoBDariaon of values of ^s-a-Jx.o) aooordlng to 
equation*) (71), (72), and (63) . T 



\^x(x,o) 


Distributed load 


(A - t 3 ) 


i \. 
t \. 


Ey « CD 


- Sy 


Oonoentxated 
load 


0 


1. 


1. 


1. 


0.016 


.994 


1. 


* WO 


0.163 


.938 


.998 


.822 


0.408 


.844 


.976 


.686 


0.817 


.705 


.888 


.648 


1.653 


.500 


.636 


.396 


3.266 


.294 


.369 


.264 


8.165 


•125 


.156 


.180 


16.330 


.064 


.078 


.065 



A 

H 

ri- 
al 

H 
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Figs. 1,2, 




(a) Symmetrical axial loading. (b) Antisymmetrical axial loading. 
Figure 1.- Sketches of panels indicating the notation used. 




(a) Completely clamped along (b) Cover sheets not attached 

edge at x = I . to support. 



Figure 2.- Sketches of oantilever box beams. 




(a) (b) 



Figure 3.- Sketches of half-sheets loaded by concentrated and distribu- 
ted forces. 
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Fig. 4 




PANEL NO. I. 

(a) 



— F 




PANEL NO. Z. 
(b) 




PANEL NO. 3. 
(C) 



A. 



Ft 

PANEL NO. 4. 

(d) 



--F 



A* 
A m 
A. 



PANEL NO. 5. 



2F 



A* 



A. 



PANEL N0.6. 
(f) 



-2F 



A. 

2A e 



iA ( 



PANEL NO. 7. 



A, 

2A. 

A. 



PANEL NO. 8. 
(h) 



2F 



A. 
2A. 



PANEL NO. 9. 

(i) 




Figure 4.- Sketches of panels snowing notation and dimensions. 




I 




Figure 11.- Stresses in cover Bheets of "box "beams with and without cut- 
out. Z — 5w. 
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figure 13.- Spanwise variation of cover- sheet effective width, for "box "beams with and without 
cut-outs . I = brt„ 
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iFigure 15,- Comparison of normal stress distributions in stiffened and unstiffenel half-sheets. 




0 * .5 1.0 1.5 x 3.0 .2.5 3.0 



v 

(a) a * 1/2. 

Figure 16.- Comparison of stress distributions predicted "by the exact procedure and 
"by the Euhn procedure for panel 8, 
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T 

Figure IV.- Comparison of stress diotrioutions predicted by the exact 
procedure, the Reissner procedure, and the Euhn procedure 
for a box "beam, v/ithout cut-out. B = lj I- 5w. 
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Fig. 18 
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Approximate solution 
Exact solution 
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1.0 



Figure 18.- Comparison of exact stresses with results obtained by re- 
taining only one term of the series in the exact solution 
for a box beam v/ith cut-out. 6 - 1J L — 5w. 
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Z 

figure 19.- Comparison of Kuhn solution with exact solution 
for a "box "beam with cut-out, 8 = 1", 1- 5w. 



